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In this work we introduce a new plateau-like inflationary model including a quadratic scalar 
potential coupled non-minimally to gravity. This potential has a dominant constant energy density 
at early times which can realize successful inflation. It also includes an infinitesimal non-zero term 
Vo responsible for explaining dark energy which causing the universe to expand accelerating at the 
late time. We show that this model predicts small tensor-to-scalar ratio of the order of r 0.01 
which is fully consistent with Planck constraints. Using the lower and upper bounds on reheating 
temperature, we provide additional constraints on the non-minimal coupling parameter £ of the 
model. We also study the preheating stage predicted by this kind of potentials using numerical 
calculations. 


I. INTRODUCTION 

In physical cosmology, cosmic inflation is the extremely rapid exponential expansion of the early universe, driven 
by a negative-pressure vacuum energy density (l]-0|. Inflation is an elegant way to resolve the horizon, flatness, 
and monopole problems in standard cosmology. Usually a scalar field with a nearly flat potential, called inflaton, is 
introduced to produce such accelerating exponential expansion. In this manner, inflaton field slowly rolls down its 
potential and then automatically produces an inflation phase. The inflationary stage is controlled by the slow roll 
parameters 
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where n = Mp 1 = (87TG) 1 / 2 , V = V{(jf) is the inflaton potential and the prime denotes derivative with respect to 
cj>. During inflation the slow roll parameters are small e, rj <C 1 and inflation ends by the breakdown of these slow 
roll conditions. For single field inflationary models one can calculate the spectral index n s and the tensor-to-scalar 
ratio r in terms of slow roll parameters as n s = 1 — 6e ± 2 rj and r = 16e, respectively. The recent Planck data 
estimate the scalar spectral index to be n s = 0.968 ±0.006 and establish an upper bound on the tensor-to-scalar ratio 
at r < 0.09 0 ]. By using Planck±BKP data, quadratic inflationary model and natural inflation are ruled out but 
models including a simple linear potential and fractional-power monomials and also models that have sub-Planckian 
held evolution are however still allowed by the data 0,0. In general, single held inflationary models with plateau-like 
shapes are more favored by Planck±WMAP data (6|. However, there are still a large number of possible models 
satisfying the Planck±WMAP constraints Q- 

It is clear that inflation has to end before nucleosynthesis for successful structure formation. Otherwise, inflation 
would leave behind a cold universe empty of matter frozen in a state of low entropy. An era after inflation called 
reheating resolves this problem. During this stage the inflaton particles coupled to other particles can decay pertur- 
batively and entropy production occurs ;2j. As the perturbative approach does not involve the coherent nature of 
inflaton field, a non-perturbative particle production known as preheating phase was then proposed &®. In this 
phase particles are produced by parametric resonance mechanism in which the inflaton held begins to oscillate near 
the minimum of its effective potential and eventually produces many elementary particles. In this sense, in the second 
stage of inflaton decaying called reheating the produced particles interact with each other and come to a state of 
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thermal equilibrium with some final temperature T re which usually corresponds to initial temperature of radiation 
epoch. Recently, it has been recognized that the reheating stage can constrain the inflationary models by calculating 
the reheating temperature as a function of spectral index n s [14M19I ] . 

In this work we introduce a new model for inflation, where a quadratic scalar potential (Eq.@) has been augmented 
with a non-minimal coupling £ to the Ricci scalar fEa. (fTIT) I which satisfies the constraints of recent CMB data and 
reheating analysis. This model has another interesting properties such as a remnant energy density appropriate to 
describe the recent accelerated expansion of the universe which has recently received attention [20l - |29j ] . We start with 
a quadratic potential coupled non-minimally to the gravity in the Jordan frame and then move to the Einstein frame 
by employing a special conformal transformation. This transformation brings the gravitational and kinetic terms in 
the Lagrangian into canonical form and also gives the potential term a plateau-like shape. Then we apply the method 
recently developed in fl4l - |T9j to put constraint on the non-minimal coupling £ of the model. In this approach, we 
calculate the reheating temperature T re and reheating e-folding number N le in terms of n s , effective equation of state 
w re and £. Using the CMB constraint on n s , the lower bound on the T re by primordial nucleosynthesis (BBN) and 
the upper bound on T re by the scale of inflation, it is possible to find a physical range for £. Indeed, the model makes 
distinct predictions for the parameters n s and rasa function of £. At the next step, using the constraint placed on £ 
we compare the n s — r prediction of the model with the recent Planck results. We show that our model is consistent 
with the Planck data well. 

In the rest of the paper, we study the preheating mechanism of the model using a four-legs interaction between 
background inflaton scalar field and a quantum scalar field through non-linear preheating mechanism. We use and 
modify the LATTICEASY program for numerical simulations [30( to understand the microphysics in our system 
during perheating era. For a four-legs interaction between the noted scalar fields, the decay of the inflaton field is not 
complete. This mechanism which was first proposed in Q, involves the idea of incomplete inflaton decay that leaves 
part of the oscillating inflaton field decoupled, behaving as dark matter. 

The paper is organized as follows: In section II, we introduce the model and discuss its motivation. In section III, 
we calculate the power spectrum of the model and check its consistency with Planck’s new data. Also, we study the 
reheating constraints on the dynamics of its potential. In Section IV we study the preheating mechanism of the model 
using a four-legs interaction term. Finally, the conclusions of the present work are summarized in section V. 


II. BUILDING A NON-MINIMAL COUPLING MODEL IN THE EINSTEIN FRAME 

The Planck data @ shows that we live in a simple universe with small spatial curvature possessing a simple dynam¬ 
ical mechanism which caused the smoothing and flattening of the universe. These properties eliminate complicated 
inflationary models with multiple fields and non-canonical ones and in exchange justify single scalar fields. 

On the other hand, we know observationally that the expansion of the universe is accelerating. The simplest ex¬ 
planation is supposing an extremely small but non-zero positive cosmological constant, A which maybe introduced 
via the vacuum state of a scalar field with non-zero constant energy Vo called dark energy. In string theory, there 
are a large number of vacua with zero or non-zero values V (< f >) in which <j> is the scalar field. Many of infinitesimal 
non-zero vacua which form a continuum lead to the creation of a large, observable universe via a self-reproducing 
inflationary mechanism. The space of all such vacua is called the landscape a There is a perfectly flat plain in the 
landscape called supermoduli-space where the vacua are super-symmetric and a quantum field potential V{4>) is zero. 
However, there are also possible non-zero values V(<p) which are stable or meta-stable local minima. Once, one moves 
off the supermoduli-space, the low energy properties of string theory which can be approximated by field theory may 
break down in some regions of the landscape and give a minimum energy for Vo which may be interpreted as the 
cosmological constant. In this way, anthropic principle enforces that we live in a special part of the universe capable 
of galaxy formation and with infinitesimal non-zero Vo- 

Therefore, combining a single scalar inflaton field potential with Vo as 

v{4>) = Uo + \ml4> 2 , ( 2 ) 

which is adopted in references like [20l l2ll |28| , would be appropriate for explaining inflation and dark energy in a 
unified scenario. During inflation, the first term of ([2]) is negligible but at the late time that cf> <C 1, it gives us 
accelerating expansion. In order for the fossil term Vo to match the dark energy density, we have to set Vo = U*a = 
Pue — 10 -47 GeU 4 [ 32 J where a is fine-tuned to be of order a ~ 10~ 112 and U* is the energy scale of inflation @ as 

V* » (1-88 x 10 16 GeV) 4 ^ . (3) 

Cosmological evidence strongly suggests that the radiation and matter dominated eras of the cosmic evolution reside 
between two exponentially expanding phases: An early phase corresponding to cosmic inflation with a nearly constant 
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Hubble parameter and a corresponding energy scale of the order of 10 64 GeH 4 and a second phase in which the 
universe has entered recently, as indicated by distant supernova observations. The recent quasi-inflationary period, 
however, is occurring at a much lower energy scale of the order of 10 _47 GeH 4 . Besides the string theory motivation 
mentioned above, the main motivation for using the present potential is that it interpolates smoothly between these 
two hierarchical energy scales. It is interesting that such a potential originates from a surprisingly simple potential 
in the Jordan frame. However, there is still a problem with flSJ) as it includes a single scalar field potentials which is 
ruled out by recent Planck data [BJ. In what follows, we try to overcome this problem by the help of the non-minimal 
coupling scenario. 

A general theory of inflation in which a single inflaton field couples non-minimally to gravity has been extensively 
studied during past years fjjl flTH . In order to develop cosmological perturbations and calculate the power spectrum, 
it is usually useful to employ a conformal transformation technique to go from the Jordan frame to the Einstein frame 
in which the calculations are simpler. In the Einstein frame, we again get a minimally coupled theory with a new 
effective potential. Fortunately, it has been confirmed that the Jordan and Einstein frames describe the same physics 
and have equivalent predictions for the CMB anisotropies [3314381 l43 - l47| . We can therefore move to the Einstein 
frame and apply the standard curvature perturbation computations. In general, one can consider the following action 
for the inflaton non-minimally coupled to gravity with non-canonical kinetic term [48l4B3j 


Sj 


d^XyJ—g 


n(4>) 

2 K 2 
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(4) 


where R is the Ricci scalar, V((f>) is the potential term and and w(</>) are functions of the inflaton field. Moving 
to the Einstein frame by the following conformal transformation 


9^ , 

leaves the action (|2]) in the following modified form 
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where we assumed ui(<j)) = Later, we show that this choice is acceptable. Based on the universal attractor 

models, the generalized non-minimal coupling to gravity is described in Joran frame with = 1 + £/(^>), where 
/(</>) can be an arbitrary function [40L l4lj. In this paper, our aim is to represent a single field plateau-like inflationary 
model describing inflation and dark energy which is consistent with recent CMB data. As explained before, the scalar 
potential in ([21) which has a minimally coupling £ = 0 to gravity is ruled out by Planck recently. On the other hand, 
the n s — r plots in [Bj favor a special subclass of single field slow-roll inflationary models with plateau-like potentials 
and with a canonical kinetic term in the framework of Einstein gravity [6 : j. Therefore, to solve the problem of @ and 
find a plateau-like potential, we use the non-minimally coupling scenario which helps us by choosing an appropriate 
function /(0) = (■f) _1 (\/l _ + k 2 £<^> 2 — 1) which gives us 

H(<£) 2 = 1 + K 2 ^ 2 . (9) 

As we do not know the exact value of inflaton mass, we define = 2 H*k 2 £. In section III, by using reheating 
analysis we put constraint on the value of £ and consequently on the value of m Now by applying (|2)l . ([9|) and recent 
definition for m^, we find F and V as 
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When inflation starts, we have F ~ 1 in the region (j) 2 1 which gives us non-minimally coupled inflaton field in 

the Einstein frame. Therefore, the former choices for Q,{4>) and change the kinetic term of © into a suitable 
canonical form. At the end of inflation, the scalar field with the plateau-like potential CCD has an oscillatory phase 
of evolution around 4> = 0 and at the late time this potential is dominated by V((f>) ~ k*a which can be interpreted 
as remnant dark energy |2ol [HI . 

Someone may think about the similarity of our non-minimally coupled model, (EH) and conformal a-attractors [55l — 
HH- The attractor nature of these theories comes from the existence of any leading pole in the kinetic term of their 
Lagrangian which lead to nearly the same predictions of cosmological observables in the limit of a large number of 
e-folds. Both classes have a free parameter and by choosing the parameter suitably one can make the model consistent 
with CMB data. However, these two classes differ except for a special case of the a-attractors with a pole in the 
kinetic term (7) . 

In the next section we fix our free parameter £ using the constraint given by the reheating stage and then compare 
the cosmological predictions of EH with Planck data. 


III. CMB CONSTRAINTS ON THE INFLATIONARY FRACTIONAL POTENTIAL 

In this section, we compare the predictions of the potential EH for n s and r with the recent CMB data. At first 
we calculate the number of e-folding N(cf>) 



( 12 ) 


where (j> e is the value of inflaton field at the end of inflation. Integrating the expression (fl2|) in the region long before 
the end of inflation, we arrive at 


m) = f M ) 4 • 

Therefore, the slow roll parameters are also given in terms of A as follows 
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Finally the n s and r parameters are obtained as 
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and hence we find the consistency relation 
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Before comparing the model with data we discuss the possible constraints of reheating stage on the parameter £ 
using the methods proposed in iT~i] ITil . The reheating era begins once the inflation ends and the comoving horizon 
starts to increase. At this time the energy density of inflaton field is dissipated into a hot plasma with temperature 
T re . The reheating can occur suddenly or it may take several e-folding times A re . After this stage, the universe 
evolves under radiation domination regime for Ard e-folding numbers until the time of equality. Taking a mode with 
comoving momentum fc, the amount of expansion during different eras are related as follow [6l| 
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FIG. 1: Plots of T re (a) and N re (b) as a function of n s in the la range from Planck 0]. In each plot, there are two allowable 
regions between blue dashed curve and red solid one which determine the upper and lower band of T re and N re for w Te = 0 
(light blue region) and uj re = 1/6 (gray region). 


where the subscripts refer to different eras, including the horizon exit (fc), reheating (re), radiation (RD), radiation- 
matter equality (eq) and the present time (0). During inflation we have H 2 = Tr 2 rA s /(2n 2 ) with the primordial scalar 
amplitude ln(10 lo A s ) = 3.089 [32]. Using the continuity equation one can write 


iVre 
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3(1 + w re ) 


In 



( 20 ) 


where w re is the reheating equation of state parameter and p e nd = (1 + d)T 4 n d is the inflaton energy density at the end 
of inflation with <5 the ratio of the kinetic energy to the potential energy. The results of the following calculations and 
their related plots are not sensetive to value of 5. Also we have p re = (n 2 / 30) g le T^ e with g le being the effective degrees 
of freedom of relativistic species upon thermalization. On the other hand, one can relate the reheating temperature 
T re to the present CMB temperature To through 


T/e = / 43 \ 1/3 a 0 a eq 

To \ll<7s,re/ CL e q &re 


( 21 ) 


where g s ,re is the effective number of light species for entropy at reheating. Combining equations (11911 and m and 
using m, we obtain the e-folding number of reheating for our model as 
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Moreover, putting the definitions of p le and p e nd into (1201) one can calculate the reheating temperature in the form 


45(1 + 6)rA s (2 £) 1 / 3 


1 + ( 2 £) 1/3 _ 
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exp 


_ ^(1 + Ur 
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(23) 


We set the values g re = g St re = 100 and plot (E3l) and (1221) in Fig. 1 (a) and (b) in terms of n s with the range within 
the la region allowed by the Planck results around n s = 0.968 ± 0.006 0- From upper and lower limits on T re we 
find that a value of u] re out of the range 0 < w re < 1/6 is disfavored from model building. Therefore, in these plots 
the results are depicted only for u> re = 0 and for w re = 1/6. If we fix the upper bound on T re to be the scale of GUT 
~ 10 16 GeV and the lower bound on T re to be the scale of BBN ~ 10 -2 GeV [§2] then the physical range for the 
non-minimal coupling parameter £ is restricted by 0.1 < £ < 0.6. Hence, the width of each curve in either Fig. 1 (a) 
and (b) corresponds to this bound on the parameter £. Our reheating results of the fractional potential model for 
the central value n s ~ 0.968 shows that in the case of uj re = 0, a number N re < 20 of e-folds and T re > 10 10 GeV is 
required, while for ui re = 1/6 a number N re < 37 of e-folds and T re > 10 2 GeV is favored. 
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FIG. 2: Constraints on r for our fractional single field potential, using Planck TT+lowP (red region), BICEP2/Keck Ar- 
ray+Planck TT+lowP (BKP) (green region) and lensing+ext (blue region) which include B-mode polarization results [j. Pink 
solid line (£ = 0.6) and yellow one (£ = 0.1) show lower and upper approximate n s — r relationship for this model assuming 
50 < N < 60. It is interesting that this model favors upper limit on ro .002 < 0.11 very well. 


In the next step, using the bound obtained on the parameter £ we can compare the predictions of the fractional 
potential model for the n s (£) and r(£) calculated in (fl6l) and 1171) with the recent Planck data. Fig. 2 shows that 
the prediction of this potential is in very good agreement with the n s — r plane favored by the Planck data Si. 
This agreement is achieved for the bounds £ = 0.1 and £ = 0.6 which were consistent with the reheating constraints. 
As one can see, the model predicts 0.004 < r < 0.012 for the gravitational waves which may be detected in future 
experiments. 


IV. PREHEATING PHASE 

As it was pointed out in section I, at the beginning of the reheating era the inflaton field begins to oscillate near 
the minimum of its effective potential and consequently approaches a sinusoidal solution with decreasing amplitude. 
Reheating was initially studied by employing the usual perturbative approach [63l [64!]. However, it was shown by 
[§!:, j§|, Hi] that the perturbative analysis does not take into account the coherent nature of the inflaton field. In general 
the decay of inflaton field can be started earlier in a preheating stage where inflaton decays non-perturbatively. In 
this phase, particles are produced through the parametric resonance mechanism @111. In this way, the coherent 
effects of the scalar field can lead to the decaying of the homogenous field much faster than the perturbative method. 
In order to study the preheating mechanism we consider a four-legs interaction of the form 

r 1 2 ±2 2 

Ant = —S <t> X , 


( 24 ) 










7 


where g is a dimensionless coupling constant and x is a quantum scalar field coupled to the classical inflaton field. In 
the external inflaton background one can decompose the quantum field \ into 

X(t, x) = J d 3 k (x* k (t)4e ik ' x + X k(t)aie-° k x ) , (25) 

where dk and a\. are the creation and annihilation operators, respectively. The Klein-Gordon equations of motion for 
the coupled fields (j) and \ in an expanding FRW universe are 


4> + 3Hj>- 4v 2 </> + V' + g 2 xU = 0 , 

a z 

(26) 

X + 3-ffx - 4 v2 X + 9 2 4> 2 X = 0 , 

(27) 


and the evolution of scale factor is given by Friedmann equations. The equations of motion of the both fields are 
oscillatory equations with effective comoving frequency w 2 = k 2 + m 2 ff a 2 . Taking m x -C during inflation, 
the effective masses of the fields will be m 2 eff = m 2 + g 2 (x 2 ) and m 2 eff = g 2 (</> 2 ), respectively. Assuming an 
inflaton field with mass oscillating with the amplitude $(f), one can usually define the dimensionless parameters 
q = (g 2 $(t) 2 ) /(4m 2 ) and Ak = 2 q + k 2 /(rn 2 a 2 ) characterizing the parametric amplification of Xk modes. For instance 
in an expanding universe with i;> 1 one can find a phase of broad parametric resonance where at each oscillation of 
the field <f>, the field x oscillates several times with sharp peaks in the density number of particles nk [§j. The limits 
on the coupling g 2 which have been discussed in mi lead to 10 3 < q < 10 5 . Indeed, we should find an optimal value 
for g large enough to produce highly efficient preheating, but small enough that the occupation numbers rik ~ 1 /g 2 
produce strong re-scattering. For gcj )o < m v there is a resonance with q <C 1 in a narrow instability band about k = m 
called narrow parametric resonance. In this regime of resonance, for each oscillation of the field </> the growing modes 
of the field \k oscillate one time. For g > 1, i.e. for oscillations with a large amplitude </> 0 , the resonance occurs for 
a broad range of k. Therefore, the broad parametric resonance is extremely more efficient |8J. 

In this section we turn to study the preheating mechanism of the potential m through the broad parametric 
resonance regime. In our case when cf> — > 0 and hence inflaton oscillates around its minimum, we can expand m in 
the form 


V(<f>) ~ V* (a + k 2 ^ 2 - k 4 ,£V + • ■ ■) , (28) 

where a is infinitesimal and 1 during reheating era. Hence, we keep only the second and third terms for our 

preheating calculations. This tree-level potential [§ 6 j has a stable minimum at <f> = 0 and the system eventually settles 
down in this ground state during time evolution. Now, by using the potential (12811 we can study the dynamics of the 
fields </> and also x after inflation through non-linear preheating mechanism. Here we modify LATTICEASY program 
[30| which solves (17B1) and (1271) and also Friedmann equation on a three-dimensional lattice using finite differencing for 
spatial derivatives and a second order staggered leapfrog algorithm for time evolution. In our choice, the comoving 
edge size of the cubic lattice is L = 20/m, number of points along each edge is N = 256, and the size of time 
step of calculation is mdt = 0.001. For the rest of this paper all simulations are performed with <j>o = 0.139 Mp, 
m<f, = 10~ 6 Mp and g 2 = 2.5 x 10 ~ 6 which is the optimal value. The £ was constrained using the reheating constraint 
in the previous section. Therefore, we set £ = 0.1 for our simulation. 

At the first step, we consider the total energy densities, e of the both fields (f> and x containing kinetic and the 
gradient energies plus their interaction term. Fig. 3 compares the contributions of these energy components of our 
potential (1751) with equivalent terms of the quadratic potential m 2 (f> 2 / 2. For the chaotic potential we have considered 
m = 10 ~ 6 Mpi. In this plot the black curves correspond to the potential (l28l) and the red ones are for the quadratic 
potential. All of the blue curves show the energy density of interaction term, g 2 (j) 2 x 2 /“2*, in the potential (1751) . As 
one can see, two curves exhibit slightly different behavior. It is significant that the contribution from inflaton field 
is dominant even after preheating, up to approximately mt ~ 175. We find that the interaction term (blue curve), 
except for a temporary time interval 80 ^ mt ^ 150, is sufficiently small that allows us to define energy density for 
each field in the following form 


(27r) 3 a 4 


drk (jJk nk 


(29) 


where is the co-moving frequency and nk is the well defined co-moving occupation number 




+ \Xk\ 2 


1 

2 


(30) 






FIG. 3: Contributions of kinetic and gradient energies of fields <j> (left plots) and \k (right plots) in comparison with quadratic 
potential. The black curves belong to our model while the red ones relate to quadratic model. Also the blue one in each plot 
represent the energy density of interaction term of our potential. 


with X k (t) = a 3/2 Xk EH- There are three stages for the dynamics that govern the preheating phase: linear 
dynamics, nonlinear dynamics and perturbative dynamics. At linear dynamical regime, through early re-scattering 
the x particles are created in the resonant band and then generation of 6<j> quanta occurs by the annihilation process 
Sxk^Xk S(f>k5(t>k■ The second stage which its evolution is under a nonlinear dynamics, starts at mt ~ 100 and 
ends at mt ~ 150. In this stage, the particle occupation number is not well-defined because the amplitude of 
interaction potential is comparable to other contributions of total energy density. At the third stage, the energy 
density contributions evolve smoothly toward the higher comoving momenta under a perturbative dynamics. 

Now we turn to the time evolution of occupation number spectra nf. and as functions of k/m. We compare the 
time evolution of these spectra of the potential (l28| with the ones of the quadratic potential in Fig. 4. In this figure 
the black curves show the nf.’ x of our potential and the red ones show the n^ x computed for the quadratic potential. 
The single green curves in each plot also show the nt’ x spectra of the potential (12811 at the final time mtf = 200. 
The spectra show the exponential growth with a resonant peak at k = fc* through parametric resonance. As we can 
see, the amplitude of the spectra for the potential (l28l) has decreased with respect to the quadratic potential. The 
spectra move down for the modes k < fc* due to the re-scattering process. The green curves show that amplitude of 
nf. decreases faster than at the end of preheating. People usually believe that preheating is followed by a reheating 
era leading to complete decay of the inflaton field. However, for a four-legs interaction, g 2 (f> 2 X 2 we applied in this 
paper, the inflaton field does not decay completely. We know that if one considers a three-legs coupling between the 
inflaton field and other quantum fields p, Qi, [ll| , the inflaton decays completely. It is also possible to consider the 
incomplete decay as a mechanism for generating dark matter [ 23 , |28|. 

V. CONCLUSION 

In this work, we introduced a new model for inflation by considering a reasonable non-minimal coupling model 
in Jordan frame and then moving to the Einstein frame through conformal transformation. The arising potential 
in Einstein frame had a plateau-like shape with a free non-minimal coupling parameter £. We studied this model 
extensively. First, we fixed the allowed range 0.1 < £ < 0.6 using the constraints on the reheating temperature for 
a favored range of 0 < oj re <1/6 in Fig. [Q Then in Fig. (2[ using the constrained £ we compared n s (£) and r(£) 
parameters predicted by this model with the recent CMB data released by Planck collaboration Q|. We found that 
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FIG. 4: Time evolution of nf. (left plot) and (right plot). The black curves are related to our potential while the red ones 
are for quadratic model. The green curves are the spectrum at the final time of our simulations. 


this model predicts small tensor-to-scalar ratio of the order of r ss 0.01 which is favored by the recent data very well. 
In the next step, we studied the out-of-equilibrium nonlinear dynamics of the interacting fields during preheating by 
the use of lattice numerical simulations. In this era, oscillations of <j> lead to parametrically resonant amplification of 
Xk and in turn the amplified modes of \k excite fluctuations of <j>. At the next stage of preheating that the interaction 
term becomes sub-dominant, we studied the evolution of the components of the total energy density and also the 
total occupation number as functions of time and compared them with equivalent terms of the quadratic model of 
inflation in Figs. [3] and [4] One can see that nf. grows exponentially faster than at the beginning of preheating due 
to rescattering and conversely its amplitude decreases faster than n £ at the end of preheating. 
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